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We present a method of subtracting the foreground contamination for the measurement of CMB 
polarization. We calculate the resultant errors on CMB polarization and temperature-polarization 
cross correlation power spectra for the high frequency instrument (HFI) aboard Planck Surveyor, and 
estimate the corresponding errors on cosmological parameters 



1 Introduction 

The upcoming satellite CMB experiment Planck surveyor offers an unprecedented opportu- 
nity to measure CMB polarization. A major hurdle in extracting the primary CMB signal 
from data, apart from noise, is galactic and extragalactic foregrounds. However, as the 
foregrounds differ from CMB in both frequency dependence and spatial distribution, one 
can hope to reduce their level in a multi- frequency CMB experiment. A multi- frequency 
Wiener filtering method to implement this scheme was developed Era and applied to cleaning 
the simulated CMB temperature map by Bouchet el al. El. They showed that the residual 
contamination after cleaning the map is much smaller than the CMB primary signal, and 
therefore the foregrounds may not be a major obstacle in the extraction of CMB temperature 
angular power spectrum. However, as the CMB polarization signal is expected to be one to 
two orders of magnitudes below the temperature signal, it is likelyjo be comparable to both 
the experimental noise and the level of foregrounds. Sethi et al. El modelled and estimated 
the level of dust polarized emission at high galactic latitudes. Dust polarization is expected 
to be the dominant contaminant for measuring CMB polarization using Planck HFI. 

In this paper, we extend the multi-frequency Wiener filtering method to include the 
polarization and temperature-polarization cross-correlation. The aim of this exercise is to 
quantify errors in estimating various power spectra and consequently the errors in cosmo- 
logical parameters. Our results are relevant for Planck HFI. 



2 Multi-frequency Wiener filtering on CMB data 



Let us denote the observed data at different frequencies as, y l v , where v indicates the fre- 
quency of the instrumental channel, and i the nature of the observed field (temperature T 
and -E-mode polarization^). J/*, for a given i, takes contributions from various galactic and 
extragalactic sources, apart from the primary CMB signal and instrumental noise. Let us 
call xp the contribution of the field j due to process p, this is the quantity we want to recover 
from the observational data y l v . We assume that there is a linear relation between them. In 
multipole space, it can be expressed as: 

yi(l,m) = A%(l,m)x J p (l,m) + bi(l,m) (1) 

where A l J p is the instrument response kernel, and b\, is the detector noise level per pixel for 
the full mission time (all repeated indices are summed over). 

The problem is now how to construct an optimal estimator x p . We chose this estimator 
to be linear in the observed data: 

4 = w;iyi (2) 

The reconstruction error for a given field and process is given by \(x p — x l p ) 2 \ = {e p ) 2 . W p l 
is chosen so as to make the variance of the reconstruction error minimal. The derivatives of 
the error with respect to the filters coefficients W p c ^ should then be zero. This condition can 
be expressed as: 

A^Wp^A^p/^x^xp^,) + Wp V i (b^b^,,) = A c w ,{x p ,x p ). (3) 

2.1 Implementation of foregrounds removal 

Eq. (|3|) is valid for the general case in which various processes, fields, and corresponding 
instrumental noises could be correlated. We consider here only uncorrelated processes and 
noises between different fields and channels. We allow for the correlation between the two 
fields T and E. With these conditions, Eq. (||) can be written as a system of four matrix 
equations which can be solved by substitution. 

Bouchet et al. u defined a quantity called the 'quality factor' to understand the quality 
of extraction of the signal corresponding to a given process. A straightforward generalisation 
of the quality factor, valid for multiple fields, can be written as: 

\XpXpt) 

where we have used Eq. (||) to write the second equality. Q 11 and Q 22 can readily be 
interpreted as the quality of the reconstruction of temperature and polarization maps, re- 
spectively. It should be noted that in the presence of cross-correlations, the quality factor 
of either field is better than the case without cross-correlations. Though the reconstruc- 
tion of temperature maps is only slightly changed by cross-correlation term, the quality of 
polarization reconstruction gets a big boost from the presence of temperature-polarization 
cross-correlation. However the meaning of term Q 21 (and Q 12 ) is not apparent. Much of 
the contribution to Q 12 comes from the term with W 11 , and therefore it is very close to 
the quality factor for the extraction of temperature and is nearly independent of the po- 
larization noise. It is not surprising as it merely tells us that to optimally reconstruct the 
cross-correlation one needs to throw out the noisy data, i.e. the polarization. However, 
the quantity of interest for us is the error in the extraction of the power spectrum of cross- 
correlation, which should not be confused with Q 12 . To get a real idea of the error bars 
of the different spectra, we must define estimators of those power spectra from the filtered 
data, and compute their covariance. 

2. 2 Unbiased estimators of power spectra 

x p (Eq. ((|)) is the data obtained after performing Wiener filtering on the multi- frequency 
maps. Our aim in this section is to use this data to write an unbiased estimator of the true 





Figure 1: Power spectra and la errors for the E-mode polarization and the temperature-polarization cross- 
correlation respectively. 



power spectrum x l p . The average power spectrum of x l p can be written as: 
This can be expanded to give: 



(5) 



(6) 



Here {i, j} stand for T and E. b l p J are generalized versions of the instrumental noise, con- 
taining contributions of the noise from different channels as well as some power leakage of 
the other processes and the other fields. We then define unbiased estimators of the spectra: 
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and compute their covariances: 
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It has been assumed that both CMB and foregrounds are Gaussian fields in the computation 
of covariances. In Fig. 1 we show the E-mode polarization and ET cross-correlation power 
spectra for sCDM model and the expected la errors on their measurement using the speci- 
fications of Planck HFI. Apart from polarized dust, an additional contribution to polarized 
foregrounds from 40% polarized synchrotron emission has been assumed. As the synchrotron 
foreground is subdominant at HFI frequencies, it doesn't affect our results. However, it is 
expected to be the major foreground for MAP and Planck low frequency instrument (LFI). 



3 Errors on cosmological parameters 

The unbiased estimators defined above and their covariances can be used to estimate the 
precision we expect on the measurement of cosmological parameters once the foregrounds 
are filtered out for any desired experiment. 



3. 1 Fisher Matrix 



The Fisher matrix is defined as an average value of the second derivatives of the logarithm of 
the Likelihood function with respect to the cosmological parameters, taken at the maximum 



Table 1: Errors on cosmological parameters. A^a is given in absolute value. 
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c 2 


h 
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n s 


«t 


T/S 


Model 


796{fiK) 2 


0.5 


0.05 


0.0 


0.05 


0.99 


0.01 


0.1 


Worst case 


3.72 % 


1.11 % 


1.93 % 


0.032 


28.3 % 


0.40 % 


538 % 


78.0 % 


Best channel 


2.5 % 


1.07 % 


1.87 % 


0.031 


2.74 % 


0.38 % 


453 % 


67.6 % 


Wiener 


1.98 % 


0.23 % 


0.34 % 


6.610- a 


19.7 % 


0.10 % 


242% 


27.0 % 



likelihood point 

d 2 c 



One can take this matrix as an estimate of the covariance matrix of the cosmological pa- 
rameters. Following Zaldarriaga et al. a we generalize the Fisher matrix approach to the 
polarized case. It can be written as: 

e=2X,Y 1 J 

where Cov(C^,C/) stands for the covariance matrix of the power spectra estimators, and 
{X, Y} £ {T, E, TE}. We will then use the covariances computed for the unbiased estimators 
defined in the previous section to take into account the foregrounds removal in the error bars 
computation. 

With the present specifications of Planck HFI, our results are shown in Table 1 for 
a variant of CDM model. The best channel case corresponds to v = 143 GHz channel 
with no foregrounds and the worst case corresponds to the same channel with foregrounds 
added as noise. As is clearly seen, significant improvement on parameter estimation can be 
achieved if the foregrounds are removed using the Wiener filtering technique discussed in 
previous sections, except for the optical depth to the last scattering surface r. An accurate 
determination of such small values of r depends on the level of CMB polarization at I < 10, 
which in our case is dominated by polarized foregrounds which cannot be easily removed. 
And therefore, a careful analysis of polarized foregrounds at large scales is necessary to 
determine such a small value of r_ The different power spectra were computed with the 
CMB Boltzmann code CMBFASTQ. 
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